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Abstract 

In this paper we consider a multi-dimensional wave equation with dynamic boundary conditions 
related to the Kelvin- Voigt damping and a delay term acting on the boundary. If the weight of the 
delay term in the feedback is less than the weight of the term without delay or if it is greater under 
an assumption between the damping factor, and the difference of the two weights, we prove the global 
existence of the solutions. Under the same assumptions, the exponential stability of the system is proved 
using an appropriate Lyapunov functional. More precisely, we show that even when the weight of the 
delay is greater than the weight of the damping in the boundary conditions, the strong damping term 
still provides exponential stability for the system. 

Keywords: Damped wave equations, boundary delay; global solutions; exponential stability; Kelvin- Voigt 
damping; dynamic boundary conditions. 



1 Introduction 



In this paper we consider the following linear damped wave equation with dynamic boundary conditions and 
a delay boundary term: 



Uu — Ait — aAut 
u(x,t) = 0, 



0. 



u tt (x,t) = 

u(x, 0) = uq(x) 
u t (x,0) = ui(x) 
u t (x,t- t) = f (x,t- t) 



du adut , , , , , 

— (x,t) + — (X,t) + HiU t (X,t) + ^ 2 u t {x,t 

av ov 



ieO, t > , 
x e T , t>0, 

x g Ti, t > o , 

x efi, 
x e Q , 

x € T l5 t G (0,r) 



(1) 



where u — u(x, i),t>0,xGf2,A denotes the Laplacian operator with respect to the x variable, fl is a 

- - d 

regular and bounded domain of R N , (N > 1), Oft = T U Ti, mes(T ) > 0, T H Ti = and — denotes 

the unit outer normal derivative, a, [i± and \ii are positive constants. Moreover, r > represents the time 
delay and «o, iti , fo are given functions belonging to suitable spaces that will be precised later. 

This type of problems arise (for example) in modelling of longitudinal vibrations in a homogeneous bar in 
which there are viscous effects. The term Aut, indicates that the stress is proportional not only to the strain, 
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but also to the strain rate. See [5]. From the mathematical point of view, these problems do not neglect 
acceleration terms on the boundary. Such type of boundary conditions arc usually called dynamic boundary 
conditions. They are not only important from the theoretical point of view but also arise in several physical 
applications. For instance in one space dimension, problem (1) can modelize the dynamic evolution of a 
viscoelastic rod that is fixed at one end and has a tip mass attached to its free end. The dynamic boundary 
conditions represents the Newton's law for the attached mass, (see [4, 1, 6] for more details). In the two 
dimension space, as showed in [24] and in the references therein, these boundary conditions arise when we 
consider the transverse motion of a flexible membrane f2 whose boundary may be affected by the vibrations 
only in a region. Also some dynamic boundary conditions as in problem (1) appear when we assume that f2 
is an exterior domain of M 3 in which homogeneous fluid is at rest except for sound waves. Each point of the 
boundary is subjected to small normal displacements into the obstacle (see [2] for more details). This type 
of dynamic boundary conditions are known as acoustic boundary conditions. 

In the absence of the delay term (i.e. fj,2 — 0) problem (1) has been investigated by many authors in 
recent years (see, e.g., [12], [13], [14], [15], [21], [22]). 

Among the early results dealing with the dynamic boundary conditions are those of Grobbelaar-Van 
Dalsen [7, 8] in which the author has made contributions to this field. 

In [7] the author introduced a model which describes the damped longitudinal vibrations of a homogeneous 
flexible horizontal rod of length L when the end x = is rigidly fixed while the other end x = L is free to 
move with an attached load. This yields to a system of two second order equations of the form 

u tt - u xx - u txx = 0, x e (0, L), t > 0, 

u(0,t) = u t (0,t) = 0, t > 0, 

< u tt (L,t) = -[u x + u tx ](L,i), t>0, (2) 

u (x, 0) = uq (x) , u t (x, 0) = vq (x) , x 6 (0, L), 
K u(L,0) = r], u t (L,0) = /i. 

By rewriting problem (2) within the framework of the abstract theories of the so-called B-evolution theory, 
an existence of a unique solution in the strong sense has been shown. An exponential decay result was also 
proved in [8] for a problem related to (2), which describe the weakly damped vibrations of an extensible 
beam. See [8] for more details. 

Subsequently, Zang and Hu [27], considered the problem 

u t t - p {u x ) xt - q {u x ) x = 0, x e (0, l) , t > o 

with 

u (0, t) = 0, p (u x ) t + q (u x ) (1, i) + kuu (1, t) = 0, t > 0. 

By using the Nakao inequality, and under appropriate conditions on p and q, they established both expo- 
nential and polynomial decay rates for the energy depending on the form of the terms p and q. 

It is clear that in the absence of the delay term and for fix = 0, problem (2) is the one dimensional model 
of (1). Similarly, and always in the absence of the delay term, Pellicer and Sola- Morales [22] considered the 
one dimensional problem of (1) as an alternative model for the classical spring-mass damper system, and 
by using the dominant eigenvalues method, they proved that their system has the classical second order 
differential equation 

miu"{t) + dW{t) + kxu{t) = 0, 

as a limit, where the parameter m\ , d\ and k\ are determined from the values of the spring-mass damper 
system. Thus, the asymptotic stability of the model has been determined as a consequence of this limit. But 
they did not obtain any rate of convergence. See also [21, 23] for related results. 

Recently, the present authors studied in [13] and [12] a more general situation of (1). They considered 
problem (1) with /i2 = 0, a nonlinear damping of the form g (ut) = \ut\ m u t instead of \i\Ut and a 
nonlinear source term / (u) — |u| p_2 w t in the right hand side of the first equation of problem (1). A 
local existence result was obtained by combining the Faedo-Galerkin method with the contraction mapping 
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theorem. Concerning the asymptotic behavior, the authors showed that the solution of such problem is 
unbounded and grows up exponentially when time goes to infinity if the initial data are large enough and 
the damping term is nonlinear. The blow up result was shown when the damping is linear. Also, we proved 
in [12] that under some restrictions on the exponents m and p, we can always find initial data for which the 
solution is global in time and decay exponentially to zero. 

The main difficulty of the problem considered is related to the non ordinary boundary conditions defined 
on Y\. Very little attention has been paid to this type of boundary conditions. We mention only a few 
particular results in the one dimensional space [16, 22, 11, 17]. 

The purpose of this paper is to study problem (1), in which a delay term acted in the dynamic boundary 
conditions. In recent years one very active area of mathematical control theory has been the investigation 
of the delay effect in the stabilization of hyperbolic systems and many authors have shown that delays can 
destabilize a system that is asymptotically stable in the absence of delays (see [10] for more details). 

In [19], Nicaise and Pignotti examined the wave equation with a linear boundary damping term with a 
delay. Namely, they looked to the following problem: 

u tt - Au = 0, ieJl,(> 0, (3) 

where f2 is a bounded domain with smooth boundary <9f2 = Tq L)Ti. On Tq, they considered the Dirichlct 
boundary conditions. While on Ti they assumed the following boundary conditions: 

du 

— (x,t)=niUt(x,t)+iJ,2Ut(x,t-T), xeTi,t>0. (4) 
av 

They proved under the assumption 

< Mi (5) 

that the solution is exponentially stable. On the contrary, if (5) does not hold, they found a sequence of 
delays for which the corresponding solution of (3) will be unstable. The main approach used in [19], is 
an observability inequality obtained with a Carleman estimate. The same results were showed if both the 
damping and the delay are acting in the domain. We also recall the result by Xu, Yung and Li [26], where 
the authors proved the same result as in [19] for the one space dimension by adopting the spectral analysis 
approach. We point out that problem (1) has been already studied by Nicaise and Pignotti in [20] for \i\ = 
and a time-varying delay. They find the same condition as the one used in this paper when fi\ = by a 
different way. However our result and our Lyapunov functional are slightly different here. (See Remark 3.2 
for more details), and we want to point out that this paper may be viewed as a continuation of the work of 
Nicaise and Pignotti [20] in which an additional damping term acts on the boundary and the study of the 
competition between these two damping terms is very interesting. 

As it has been proved by Datko [9, Example 3.5], systems of the form 

f w tt - w xx - aw xxt = 0, x G (0, 1), t > 0, 

(6) 

L to(0,t) = 0, w x (l,t) = -kw t (l,t-T), *>0, 

where o, k and r are positive constants become unstable for an arbitrarily small values of r and any values 
of a and k. In (6) and even in the presence of the strong damping —aw xx t, without any other damping, the 
overall structure can be unstable. This was one of the main motivations for considering problem (1). (Of 
course the structure of problem (1) and (6) are different due to the nature of the boundary conditions in 
each problem). 

Subsequently, Datko et al [10] treated the following one dimensional problem: 

Utt(x, t) — u xx (x, t) + 2au t {x, t) + a u{x, t) = 0, < x < 1, t > 0, 
u(0,*) = 0, t>0, (7) 

u x (l,t) = -ku t (l,t-T), t>0, 
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which models the vibrations of a string clamped at one end and free at the other end, where u(x,t) is the 
displacement of the string. Also, the string is controlled by a boundary control force (with a delay) at the 
free end. They showed that, if the positive constants a and k satisfy 

e 2a + l 

then the delayed feedback system (7) is stable for all sufficiently small delays. On the other hand if 

, e 2a + 1 i 

then there exists a dense open set D in (0, oo) such that for each r £ D, system (7) admits exponentially 
unstable solutions. 

As a consequence of what we have said before, two main questions naturally arise here: 

• Is it possible for the damping term —Au t to stabilize system (1) when the weight of the delay is greater 
than the weight of the boundary damping (i.e. when [ii > /Hi)? 

• Does the particular structure of the problem prevents the instability result obtained in [9] for problem 
(6)? 

One of the main purpose of this paper is to give positive answers to the above two questions. More 
precisely, we study the asymptotic behavior (as t — > oo) and related decay rates for the corresponding- 
solutions of system (1) where the question to be addressed here is whether the delay term [i2Ut(x, t — r) can 
destroy the stability of the system, which is exponentially stable in the absence of that delay [12]. As we 
shall see below, the presence of the strong damping term aAu f in (1) plays a decisive role in the stability 
of the whole system if (5) does not hold. Thanks to the energy method, we built appropriate Lyapunov 
functionals lead to stability results. 

The paper is organized as follows: in the next section, we prove the global existence of the solutions by 
using the Lumer-Phillips' theorem in the same way as in [19]. In section 3, we show that if the weight of the 
delay is less than the weight of the damping, then the energy defined by (38) decays exponentially to zero. 
We also prove that even if the weight of the delay is greater than the weight of the damping, the solution 
still decays to zero exponentially provided that the damping parameter a satisfies an appropriate condition. 
Let us mention that without the damping factor a, Nicaise and Pignotti [19] proved the instability of the 
null stationary solution in the case [ii > fj,±, whereas we will show that if [12 > , by adding a condition of 
the form a > (/12 — Hi)B 2 (with B a constant defined later), we are able to prove the stability of the null 
stationary state thanks to a suitable choice of a Lyapunov function. 

2 Well-posedness of Problem (1). 

In this section we will first transform the delay boundary conditions by adding a new unknown. Then as 
in [19], we will use the Lumer-Phillips' theorem to prove the existence and uniqueness of the solution of 
problem (1). 



2.1 Setup and notations 

We denote 



We set 71 the trace operator from H^ o (fl) on L 2 (ri) and # 1/2 (ri) = 71 (i?^ (fi)) . We denote by B the 
norm of 71 namely: 

Vu 6 Hk (SI) , ||tt|| 2 ,r 1 < B||Vtt|| 2 ■ 
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We recall that H^ 2 (Ti) is dense in L 2 (Ti) (see [18]). 

We denote E(A,L 2 (Q)) — {u G such that Am G L 2 (fl)} and recall that for a function it £ 

E(A,L 2 (fl)) , e £T 1/2 (Ti) and the next Green's formula is valid (see [18]): 

Vu(x)Vv(x)dx = ( -Au(x)v(x)dx + l^-;v) , Vv G i?p (fi), (8) 

where (.; .) r means the duality pairing between iJ _1//2 (ri) and i? 1 / 2 ^). 

By (., .) we denote the scalar product in L 2 (f2) i.e. (it,u) = / u(x)v{x)dx. Also we mean by ||.|| g the 

Jn 

L q (il) norm for 1 < q < oo, and by ||.|L,ri the ^ 9 (ri) norm. 

Throughout the paper, we use the standard notations as in the book [3] for example. 

In order to prove the local existence of the solution of problem (1), we consider the following two cases : 

case 1: /12 < fJ>i- We may define a positive real number £ such that: 

< f < r (2fl! - fj, 2 ) ■ (9) 
case 2: /12 > Mi- We will suppose that the damping parameter a verifies: 

a > ( M2 -Mi)5 2 • (10) 
In this case, we may define a positive real number £ satisfying the two inequalities: 



a > ['^ + i-- Ml ]B^>0. (12) 



£ > t M2 , (11) 
2 2r 

2.2 Semigroup formulation of the problem 

In this section, we prove the global existence and the uniqueness of the solution of problem (1). To overcome 
the problem of the boundary delay, we introduce, as in [19], the new variable: 

z (x, p, t) — u t (x, t — rp) , x G Ti, p G (0, 1) , t > 0. (13) 

Then, we have 

rz t (x, p, t) + z p (x, p, t) = 0, in Ti x (0, 1) x (0, +00) . (14) 
Therefore, problem (1) is equivalent to: 

utt — Au — aAu t — 0, x G fi, t > , 

rz t (x, p, t) + Zp(x, p, t) = 0, x G Ti, p G (0, 1) , t > , 

m(x, i) = 0, x G To, £ > , 

. . / du dut / \ , . , A 

u t t(x,t) = - I — (x,t) + a — (x,t) + piu t (x,t) + p 2 z{x, l,t) I , xeTx, t>0, ( 15 ) 

z(x, 0, t) — u t (x, t), x G Tx, t > , 

it(x, 0) = wo(ir), a; G , 

«t(:c, 0) = ui(x), x G 51 , 

p, 0) = / (a;, -rp), x G r 1; p G (0, 1) . 

The first natural question is the existence of solutions of the problem (15). In this section we will give a 
sufficient condition that guarantees the well-posedness of the problem. 
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For this purpose, as in [19], we will use a semigroup formulation of the initial-boundary value problem 
(15). If we denote V := (u, Ut, 71 (ut), z) , we define the energy space: 

je = ffr ( fi ) x L 2 ( n ) x L2 ( r i) x ^ 2 ( r i x (°> !))■ 
Clearly, is a Hilbert space with respect to the inner product 

(Vi,V2)jg>= / Vui.Vu2<ia;+ / v \V2dx 4- / w\W2da + £ / / z\Z2dpda 
Jn Jn Jti Jt! Jo 

for 14 = (uij Wi, »i, zi) T , V2 = (u2, v 2 , W2, Z2) T and £ is defined by (9) or (11). 

Therefore, if Vq 6 ^ and V G Jff, the problem (15) is formally equivalent to the following abstract evolution 
equation in the Hilbert space Jff: 

j V'(t) = sx?V{t), t > 0, 
I V(0) = V , 



(16) 



where ' denotes the derivative with respect to time t, Vq := (u 0} u\, 71 (iii), /o(-, — -t)) and the operator 
is defined by: 

( v \ 

Au + aAv 



( « \ 

v * y 



du dv 

/ZlU - /X 2 2 (., 1) 

1 



V 



/ 



The domain of is the set of V = (it, v, w, z) T such that: 

(u,v,w,z) T e fl* D (n) x Bj o (n) x L 2 (rx) x i 2 (ru^ci)) 

u + ave E(A, L\n)) , 8 ( U + H G L 2 ^) 

01/ 

w = 7x(u) = 0) on Ti 

The well-posedness of problem (15) is ensured by: 



(17) 

(18) 
(19) 



Theorem 2.1. Suppose that P2 > p\ and a > (^2 — Hi)B 2 or P2 < p>\- Let Vo G J^, then there exists a 
unique solution V £ C (R+; Jif) of problem (16). Moreover, if Vo G 3>(g/), then 

V EC (K+; 9 GO) n C 1 (R+; JT) . 

Proof. To prove Theorem 2.1, we use the Lumer-Phillips' theorem. For this purpose, we show firstly that 
the operator srf is dissipative. Indeed, let V — (u, v, w, z) T G & {^)- We have 



\Ju.\Jvdx 



du 



v (Au + aAv) dx 
Ov 



w ( - a— - fiiv - i± 2 z (cr, 1) J do" - 
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r, Jo 



zzpdpda. 



d(u -\- otv) 

But since 11 + au e E(A,L 2 (Vl)) and G i 2 ^), we may apply Green's formula (8) where the 

duality pairing (.; .) r is simply the i 2 (Fi) inner product (because w = 71 (v) G L 2 (Fi)) and obtain: 



(s^V, V)jg, = —fix / w 2 d<j — fi2 / z (cr, 1) wda 



Vi 



IVurdx- - 
r 



ri Jo 



z„zdpdx. 



(20) 



G 



At this point, we have to distinguish the following two cases: 

Case 1: We suppose that [12 < Ml- Let us choose then £ that satisfies inequality (9). Using Young's 
inequality, (20) leads to 

(siV,V)^ + a I \Vv\ 2 dx+(^-^-~^] f w 2 da 



2r 2 r 

2r 2 )/ ri ^l)^<0. 
Consequently, by using (9), we deduce that 

{■zfV,V)^<0. (21) 

Case 2: We suppose that \ii > p\ and a > (/12 — Pi)B 2 . Let us choose then £ that satisfies the two 
inequalities (11) and (12). Using Young's inequality and the definition of the constant B, we can again prove 
that the inequality (21) holds. This means that in both cases si is dissipative. 

Now we show that A J — si is surjective for all A > 0. 

For F = (/1, / 2 , f 3 , U) T 6 Jf, let V = (u, v, w, z) T e 9 {si) solution of 

{XI — si) V = F, 

which is: 



Xu-v = f u (22) 
Xv-A(u + av) = / 2 , (23) 



d{u + av) 



Xw + ^-g^ — + Mi" + A»2«(., 1) = /a, (24) 

Az+iz p = / 4 . (25) 
r 

To find U = (u, w, id, z) T G S> {si) solution of the system (22), (23), (24) and (25), we proceed as in [19]. 
Suppose u is determined with the appropriate regularity. Then from (22), we get: 

v = Xu-fi . (26) 

Therefore, from the compatibility condition on T%, (19), we determine z{.,0) by: 

z{x,0) = v(x) = Xu{x) - fx{x), for 16^. (27) 

Thus, from (25), z is the solution of the linear Cauchy problem: 



z p = r(fi{x) - Xz{x,p)\ for x G T 1 , p G (0, 1), 
z{x, 0) = Xu{x) — fi{x). 



(28) 



z{x, p) = Xu{x)er xpT - fie- xpT + re" Apr / f 4 {x, a)e XaT da for x G T 1 , p G (0, 1). (29) 



The solution of the Cauchy problem (28) is given by: 

JO 

So, we have at the point p = 1, 

-At 



z{x,l) = Xu{x)e- XT +zi{x), forxGTi (30) 
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with 



[ f 4 (x,a)e x,yT d(j, for a; G IV 
Jo 



Since /i G #r (fi) and / 4 £ i 2 (ri) x L 2 (0, 1), then zi G £ 2 (Ti). 

Consequently, knowing it, we may deduce v by (26), z by (29) and using (30), we deduce w — 71(f) by 
(24). 

In view of equations (23) and (24), we set, as in [20], u = u + av. Then, from (26), we have 

v = Xu — fx = X(u — av) — fi . 
Since A > and a > , 1 + Xa ^ 0; thus we have: 

A _ /1 , , 

v = ~« r - u _ ; r - - (31) 

1 + Xa 1 + Aa v ' 

But since u = u — av, we have: 

' - " f, (32) 



1 + Xa 1 + Xa 
From equations (23) and (24), u must satisfy: 

-r-j-T—u - Au = / 2 + 7-77— /l) infi (33) 
1 + Aa 1 + Aa 

with the boundary conditions 

u = 0, on T (34) 

— = / 3 - A7i(u) -/Xi7i(«) -^(.,1), on r\ (35) 
01/ 

the last equation at least formally since we don't have yet found the regularity of u. Replacing u by its 
expression (32) and inserting it in equation (30), we get: 



z(x, 1) = - A x u(x)e Xt + A " h{x)e Xt + Zl (x), for x G IV 
1 + Xa 1 + Xa 

Using the preceding expression of z(., 1) and the expression of v given by (31), we have 



m X(fi 2 e- XT + {X + m) 

t- = - — rn + for 2:6111 ( 36 ) 

w 1 + Aa 

with 

s , , x . (A + Mi) -/i 2 Aae~ Ar 

7W = 73 W H r— : JiW - M2^i(x), for x G ii . 

1 + Aa 

From the regularity of /1 , / 2 , zi, we get / 6 i 2 (ri). 

The variational formulation of problem (33), (34), (36) is to find u G (fi) such that: 

\a /■ A(^ 2 e- AT + (A + /i 1 ) 

777/ 1 I V7T7V7, . J„ I / > / ST 



hw + VuVwdx + / u(a)u}(a)da, (37) 



q 1 + Aa 7ri 1 + Aa 

72 + Y~j~^/i^ wcte + J /(oOw(cr)do-, 

for any ui G (fi). Since A > , /ii > , // 2 > 0, the left hand side of (37) defines a coercive bilinear form 
on Hh (fi). Thus by applying the Lax-Milgram theorem, there exists a unique u G Hf. (fi) solution of (37). 



Now, choosing oj £ < ^ c °°, u is a solution of (33) in the sense of distribution and therefore u £ E(A, L 2 (Q)). 
Thus using the Green's formula (8) in (37) and exploiting the equation (33) on 17, we obtain finally: 



L 



\U 2 e- XT + (A + Mi)_ i m 
— <a)u{a)da + ( -; U 



f(a)uj{a)da Vlo e H^ (Cl) 



So 5 £ E(A, L 2 (n)) verifies (36) and by equation (32) and (31) we recover u and v and thus by (29), we 
obtain z and finally setting w — 71 (w), we have found V = (u, v, w, z) T g 3> solution of (Id — &/) V = F. 
Thus, the proof of Theroem 2.1, follows from the Lumer-Phillips' theorem. □ 



3 Asymptotic behavior 

3.1 Exponential stability for ^2 < Hi 

In this subsection, we show that under the assumption ^2 < the solution of problem (15) decays to the 
null steady state with an exponential decay rate. For this goal, we use the energy method combined with 
the choice of a suitable Lyapunov functional. 

For a positive constant £ satisfying the strict inequality (9), (i.e. < instead of <) we define the functional 
energy of the solution of problem (15) as 



E(t) = E(t,z,u) 



Iv^ii^ + ii^wil' + II^W" 2 



2,ri 



r x Jo 



z 2 (<7, p, t) dp da 



T t JO 



z 2 (er, p, t) dp da, 



(38) 



where 



E 1 (t) = \\S7u(t)\\l + \\u t (tnl + \\u t (t)\\lr 1 - 
Let us first remark that this energy is greater than the usual one of the solution of problem (1), namely 
Ex(t). 

Now, we prove that the above energy E (t) is a decreasing function along the trajectories. More precisely, 
we have the following result: 

Lemma 3.1. Assume that p\ > then the energy defined by (38) is a non-increasing positive function 
and there exists a positive constant C such that for (it, z) solution of (15), and for any t > 0, we have: 



dE (<) 
dt 



< -C 



u t (o~,t) da + / z (a,l,t)da 



a I \Vut(x 7 t)\ dx 



(39) 



Proof. We multiply the first equation in (15) by ut and perform integration by parts to get: 



ld_ 

2di 



\\vu(t)f 2 + \\u t (t)\\l + \\ut(ty 2 



+ a\\Vu t (t)\ 



||ut(£)|| 2 r + P2 / u t (a,t)u t (a,t - r)da = . 
We multiply the third equation in (15) by £z, integrate the result over Ti x (0, 1), we obtain: 

c /• r , „ , , € r f 1 o 



(40) 



• / / z p z(a, p,t) dp da 
Jr x Jo 



2Tj ri io dp 



-z (a, p, t) dp da 



|- f (z 2 (a,l,t)-z 2 (a,0,t))da 
lT Jri 



(41) 
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Using the definition (13) of z in the equality (40) and using the same technique as in the first step of the 
proof of Theorem 2.1, where we proved that is dissipative, inequality (39) holds. □ 

The asymptotic stability result reads as follows: 

Theorem 3.1. Assume that p2 < p.\. Then there exist two positive constants C and 7 independent oft 
such that for (u,z) solution of problem (15), we have: 

E{t) < Ce-t\ Vf > . (42) 

Proof. The proof of Theorem 3.1 relies on the construction of a Lyapunov functional. 
For a small positive constant e to be chosen later, we define: 

L(t) = E(t) + e u(x,t)ut(x,t) dx + e / u(a,t)u t (a,t) da 
J n Jr ± 

+ — / \Vu(x,t)\ dx (43) 
2 Jn 



+ e£ [ [ e- 2Tp z 2 (a,p,t) dp da. 
Jr, Jo 



Let us say that the introduction of the last term in the Lyapunov functional L is inspired by the work of 
Nicaise and Pignotti [20]. 

It is straightforward to see that for e > 0, L(t) and E{t) are equivalent in the sense that there exist two 
positive constants f3\ and fa > depending on e such that for all t > 

faE{t) < L(t) < faE(t) . (44) 

By taking the time derivative of the function L defined by (43), using the equations in problem (15), several 
integration by parts, and exploiting (39), we get: 



dL(t) < _ c 



uf(a,t) da + / z (a, l,t) da 
r, Jri 



dt 

- a\\Vu t \\l - s\\\7u\\l + e\\u t \\l + e\\u t \\l ri 

— e/ii / u t (a,t)u(a,t) da — Efi2 / z(a,l,t)u(a,t)da (45) 



+ £ tj t (J r I e- 2Tp z 2 (a,p,t)dpda^ 
By using the second equation in (15), the last term in (45) can be treated as follows: 



-2rp 



z(a, p, t)z p {a, p, t) dp da 



= -— I ^ e- 2Tp -^-z 2 (a,p 7 t)dpda. 
t Jti Jo d P 

Then, by using an integration by parts and the definition of z, the above formula can be rewritten as: 



*^{L Si e_2TPz2((T '^ ) dp da 



—e- 2r [ z 2 (l,p,t) dpda+— I u 2 (a,t) da (46) 
T Jri t 



-2e£ [ f e- 2Tp z 2 (a,p,t) dp da 
Jrx Jo 
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Applying Young's inequality, and the trace inequality, we obtain, for any 8 > 0: 



Utuda 



< <y||u||2,r l + 4£lkll2,r 1 

< 5B 2 \\\7u\\l + ±- s \\u t \\l ri . 



Similarly, we have 



z{a, 1, t)u{a, t) da 



-Is J z2 ( a ^^) da + 6B2 W S7u \ 



(47) 



(48) 



Inserting (46), (47) and (48) into (45) and using Poincare's inequality for u t , in which we denote C(tt) the 
Poincare's constant, namely : 

\/w e flj (Q) , ||w|| 2 < C(n)||V«;||a 



we have: 



dL(t) 
dt 



< 



C-e 1 



£ Mi 
t AS 



\ u th,ri 



/ z 2 (a, 1, t) da 



a - eC(n) 2 ) \\Vu t \\l -e(l-B 2 5 ( Ml + M2 ) ||Vu||2 



-2e£ / / e~ 2Tp z 2 (a,p,t) dp da. 
Jr! Jo 



(49) 



We choose now S small enough in (49) such that 

5< 



B 2 (m+fi 2 ) ' 



Once S is fixed, using once again Poincare's inequality in (49), we may pick e small enough to obtain the 
existence of r) > 0, such that: 

dUt) 



dt 

On the other hand, by virtue of (44), setting 7 = — r)e/p2, the last inequality becomes: 

dL(t) 



dt 



< -jL(t) , Vt > 



(50) 



(51) 



Hence, integrating the previous differential inequality (51) between and t, we get 

L(t) < C*e~ 7t , Vt > 0, 

for some positive constant C* . 

Consequently, by using (44) once again, we conclude that it exists C > such that: 

E{t) < Ce" 7 * , Vt > . 



This completes the proof of Theorem 3.1 



□ 
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3.2 Exponential stability for ^ > £*i an d a > {1x2 — ^i)B 2 

As, we have said in the Introduction, and it is clearly observed in Theorem 3.1, that the strong internal 
damping compensates the destabilizing effect of the delay in the boundary condition. 

In this section, we assume that 1x2 > Mi and a. > (/12 — ii\)B 2 . As we will see, we cannot directly perform 
the same proof as for the case where /12 < >Ui , since the boundary delay term fj,2Ut(x,t — r) is greater than 
the normal one /iiUt(x, t), i.e. (/12 > l^i). So we have to control this term by the damping term aA« ( in the 
equation. 

Remark 3.1. In the case 2, namely [12 > /ii the condition a > (p2 — ^i)B 2 permits us to find £ satisfying 
(11)-(12). This choice of £ is essential in the proofs of Lemma 3.2 and Theorem 3.2 below. 

Lemma 3.2. Assume that /i 2 > [i\ and a > (^ 2 — ll\)B 2 . For any £ satisfying (11)-(12), the energy defined 
by (38) is a non-increasing positive function and there exists a positive constant k such that for (u, z) solution 
of (15), and for any t > 0, we have: 



dEit) < k 



z 2 (a, l,t) da + / \Vu t (x,t)\ 2 dx 



dt 

Proof. Let us first recall from (40) and (41) the following identity 
dE(t) 



(52) 



dt 



a J \Vu t \ 2 dx - ^/xi - J u 2 t {a,t) da 
— [ z 2 (a, 1, t) da — /i 2 [ u t (a,t)z(a,l,t) da 

2r JTi JT! 



(53) 



Now, using Young's inequality, then (53) takes the form: 



dE(t) < -af Q |V Ut | 2 dx - fa - |- - f ) jf u 2 (a,t) da 



dt 



JL _ t± 
2t 2 



z 2 {a,\,t) da. (54) 
ri 



Since, 



£ ^ 2 ^ n 



then, using the trace inequality, we obtain: 



z 2 (a, 1, t) da. 



dt ~ V V 2t 2 

JL _ ^£ 
> 2 / Jri 

Using the two inequalities (11)-(12) that £ satisfy , we may find n > such that the inequality (52) holds. □ 

We can now state that under the same assumption as in Lemma 3.2, the system (15) is also exponentially 
stable. The second stability result reads as follows: 

Theorem 3.2. Assume that /i 2 > ll\ and a > (/i 2 — ^i)B 2 . For any £ satisfying (11)-(12), there exist two 
positive constants C andj independent oft such that for (u,z) solution of problem (15), we have: 

E(t) <Ce-^, Vt > . (55) 
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Proof of Theorem 3.2. We use the same Lyapunov function as in the previous section, namely, for a small 
positive constant e to be chosen later, we define: 



L(t) = E(t) 



u(x,t)ut(x,t) dx + e I u(a,t)ut(a,t) da 
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ea 

Y 



|Vu(x, f)| dx 



r x Jo 



-2rp 2 



z (a, p, t) dp da. 



By taking the time derivative of the function L, using the equations in problem (15), several integration 
by parts, and exploiting (52), we get: 



dL(t) 
dt 



< 



z 2 (a,l,t)da + / |Vi/t(x, f)| dx 



a\\Vu t \\l ~ e||V«||l + e\\u t \\l + e\\u t \\ 



2Ji 



Ut(a,t)u(a,t) da — ep2 / z(a,l,t)u(a,t)da 



(56) 



dt 



ri Jo 



-2rp„2 



z (a, p, t) dp da 



Let us remark that equality (46), and inequalities (47) and (48) used in the the proof of Theorem 3.1 are 
still valid. Inserting (46), (47) and (48) into (56) and using Poincare's inequality, we have: 



dL(t) 
dt 



< 



/ z 2 (a,l : t)da 



K-e B 



t 45 

(i-c(n) 2 s(p 1 + p 2 ))\\wu\\ 



|V«f 



(57) 



ri Jo 



-2t P 2 



z (a, p, t) dp da . 



The remaining part of the proof is similar to the one of the proof of Theorem 3.1: by choosing firstly S and 
then e, we may find 7 > independent of t such that: 



dL(t) 
dt 



< -jL(t) , yt > 



This inequality permits us to conclude the proof of Theorem 3.2. 



□ 



Remark 3.2. After our work had been submitted, we noticed that a similar problem has been already studied 
by Nicaise and Pignotti [20]. However our problem is slightly different from the one considered in [20]: 

• First, we have the extra term p±Ut on the boundary conditions, which makes the stability analysis 
independent of the strong damping aAu t , for p\ > pi- This is not the case in [20]. See the assumption 
(2.56) in [20]. Moreover, the paper [20] does not cover the case a — at all. 

• Secondly, our Lyapunov functional is different from the one used in [20]. We choose to use this one to 
have a strong control of the boundary delay term. 

Finally, let us remark that our assumption (/i 2 — Pi)B 2 < a is exactly the same than the one obtained in 
[20, Condition (2.56)], when p\ = 0. So this work can be viewed as a continuation of the works of Nicaise 
and Pignotti [20] 
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